Understanding Anomalous Transport in Intermittent Maps: 
Prom Continuous Time Random Walks to Fractals 
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We show that the generalized diffusion coefficient of a subdiffusive intermittent map is a fractal 
function of control parameters. A modified continuous time random walk theory yields its coarse 
functional form and correctly describes a dynamical phase transition from normal to anomalous 
diffusion marked by strong suppression of diffusion. Similarly, the probability density of moving 
particles is governed by a time-fractional diffusion equation on coarse scales while exhibiting a 
specific fine structure. Approximations beyond stochastic theory are derived from a generalized 
Taylor-Green-Kubo formula. 
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The notion of anomalous diffusion derives from the 
fact that the mean squared displacement (MSD) must 
not follow the law of normal diffusion, < x 2 >~ t 1 with 
7 = 1, but may be supcrdiffusive with 7 > 1 or subdiffu- 
sive with 7 < 1. Such anomalous dynamics has been in- 
vestigated theoretically and observed experimentally not 
only in amorphous semiconductors, surface diffusion, tur- 
bulence, polymers and plasmas but also in chemical, bi- 
ological and economical problems [1, 2]. For this vari- 
ety of applications it is desirable to have a class of test 
systems at hand which are yet easy to handle. Here pe- 
riodically continued deterministic maps provide an im- 
portant basis for analytical and numerical investigations 
[3, 4, 5, 6, 7, 8]. Both the sub- and the supcrdiffu- 
sive dynamics of such models was successfully described 
by continuous time random walk (CTRW) approaches 
[3, 5, 6, 8]. This stochastic theory yields non-Gaussian 
probability density functions (PDF) for anomalous dif- 
fusive processes exhibiting stretched exponential decay 
in case of sub- and Levy power law decay in case of su- 
perdiffusion [5]. 

On the other hand, a microscopic theory of anomalous 
deterministic transport explaining the origin of this be- 
havior in terms of thc theory of dynamical systems is just 
beginning to evolve [7]. Surprises along these lines were 
already encountered in normal diffusive maps, where the 
diffusion coefficient was found to be a fractal function of 
control parameters [9, 10]. This behavior was also re- 
ported for more complex systems like the climbing sine 
map [11], bouncing ball billiards [12] and coupled Joseph- 
son junctions [13]. The fractality can be understood as a 
signature of long-range dynamical correlations that, due 
to topological instabilities, change in a complicated man- 
ner under parameter variation [14]. However, no such 
fractal structures were yet identified in anomalous dy- 
namics as modeled by the maps of Refs. [3, 4, 5, 6, 7, 8]. 

In this Letter we focus on a subdiffusive map whose 
functional form on the unit interval was introduced by 



Pomeau and Manneville for describing intcrmittency [15], 



x n+1 = M a , z (x n ) = x n + ax z n , <x n < 



1 



(1) 
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< x n < 1 the map is defined according to 
M a)Z (— x) = —M atZ (x). The translation M Q>z (a; + m) = 
M atZ (x) +771, m £ Z, completes the definition on the real 
line, see the inset of Fig. 1 (a). The degree of nonlinear- 
ity z > 1 and a > 1 in Eq. (1) hold for the two control 
parameters. For 1 < z < 2 this model leads to nor- 
mal diffusion while for z > 2 its behavior is anomalous 
[3, 5]. The anomalous regime results from the existence 
of marginal fixed points located at all integer values of 
x. Thus, a typical trajectory of the map consists of long 
laminar phases interrupted by chaotic bursts. In contrast 
to Refs. [3, 5], which focused on the time dependence of 
the MSD for the particular parameter value of a = 2 Z , 
here we study the behavior of the generalized diffusion 
coefficient (GDC) [2] 
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(x 2 ) 

lim 

n — too fCi 
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(< ... > denotes an ensemble average) for this map under 
variation of both control parameters. We first focus on 
numerical simulations [16] of the map Eq. (1) and on the 
interpretation of the results within the CTRW approach. 
CTRW theory models diffusion processes by sequences of 
jumps interrupted by periods of waiting. Let the PDFs 
of waiting times and of jump lengths be defined by 4>{t) 
and A(x), respectively. Choosing these functions appro- 
priately, CTRW theory predicts that 7 = 1 for 1 < z < 2 
and 7 = l/(z — 1) for 2 < z irrespective of the parameter 
a of the map [3, 5]. Indeed, for all a we find excellent 
agreement between these solutions and the values for 7 
obtained from simulations. Consequently, the analytical 
expressions for 7 are used throughout our work. 

Simulation results for K as a function of a for fixed z 
are presented in Fig. 1. Magnifications of Fig. 1 (a) shown 
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in parts (b) and (c) reveal self similar-like irregularities 
indicating a fractal parameter dependence of K. Note 
particularly the structure marked by triangles, which is 
repeated on finer and finer scales. The parameter values 
for these symbols correspond to specific series of Markov 
partitions [9, 10]. In Fig. 2 K is depicted as a function 
of z for different values of a, again displaying highly non- 
monotonic parameter dependences. 

In a first step for explaining these curves we apply stan- 
dard CTRW theory [3, 5] by modifying this approach at 
three points: Firstly, the waiting time PDF must be cal- 
culated according to the grid of elementary cells indicated 
in Fig. 1 [10, 17] yielding 

4>{t) =a(l + a(z- l^r^ 1 ■ (3) 
Secondly, as a jump probability we use 

\(x) = ^S(\x\-l) + (l-p)5(x). (4) 

The second term accounts for the fact that particles may 
stay in a box with a probability of (1 — p). This quantity 
is roughly determined by the size of the escape region p = 
(1 — 2x c ) with x c as a solution of the equation x c + ax z c = 
1. Thirdly, we introduce two definitions of a typical jump 
length k,ie {1,2}. Thus, 

h = {\M atZ {x)-x\} (5) 

corresponds to the actual mean displacement while 
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FIG. 1: The generalized diffusion coefficient K, Eq. (2), as 
a function of a for z = 3. The curve in (a) consists of 1200 
points, the dashed-dotted line displays the CTRW result K\, 
Eqs. (5), (7). (b) (600 points) and (c) (200 points) show mag- 
nifications of (a) close to the onset of diffusion. The dotted 
line in (b) is the CTRW approximation K2, Eqs. (6), (7). 
The dashed line represents the first term of the TGK formula 
Eq. (8). The triangles mark a specific structure appearing on 
finer and finer scales. The inset depicts the model Eq. (1). All 
quantities here and in the following figures are dimensionless. 




FIG. 2: K as a function of z for a = 3.14,4,5 consisting of 
200, 250 and 300 points (bottom to top). Bold lines represent 
simulation results for n = 10 4 iterations, the dashed line is 
the CTRW solution K 2 , Eqs. (6), (7), for n -> 00. The thin 
line displays the first term of the TGK formula Eq. (8). At 
z — 2 there is a transition from normal to anomalous dif- 
fusion. The inset describes the time-dependence of < x 2 > 
around this value with a — 3.14. Here dashed lines correspond 
to z = 1.5, 1.9, 1.95 (bottom to top), dashed-dotted lines to 
2.5, 2.1, 2.05 (top to bottom) while the bold line depicts z = 2. 
The dotted line is proportional to ln(n). 



i 2 = {|[M 0l ,(ar)]|} (6) 

gives the coarse-grained displacement in units of elemen- 
tary cells, as it is often assumed in CTRW approaches. In 
these definitions {. . .} denotes both a time and ensem- 
ble average over particles leaving a box. The modified 
CTRW approximation for the GDC resulting from these 
three changes reads 

R = fp|?oT sin(7r 7 )/7r 7 1 +T, < 7 < 1 (?) 
1 \plja(l - I/7), 1 < 7 < co . U 

Fig. 1 (a) shows that K\ well describes the coarse func- 
tional form of K for large a. K2 is depicted in Fig. 1 
(b) by the dotted line and is asymptotically exact in the 
limit of a — ► 0. Hence, the GDC exhibits a dynami- 
cal crossover analogous to the one found for normal dif- 
fusion [10, 11, 17]. Let us now focus on the GDC at 
a = 12, 20, 28, which corresponds to integer values of 
the height h = [M a , z (l/2)] of the map. Simulations re- 
produce, within numerical accuracy, the results for K2 by 
indicating that K is discontinuous at these parameter val- 
ues. Due to the self similar-like structure of the GDC we 
arrive at the conjecture that K exhibits infinitely many 
discontinuities on fine scales as a function of a, in contrast 
to the continuity of the CTRW solution. 

We now switch back to the z-dependence of K, where 
CTRW theory predicts a dynamical phase transition from 
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normal to anomalous diffusion at z = 2 [3, 5]. At the 
transition point one has < x 2 >^ n/log(n). Accord- 
ing to Eq. (2), Eq. (7) yields a continuous transition 
in z, see K 2 in Fig. 2. However, this coarse functional 
form is obscured not only by fractal irregularities but also 
by an ultra-slow convergence of the MSD in our simula- 
tions. Using CTRW theory we verified that around z = 2 
the MSD is determined by a series of logarithmic correc- 
tions in time, where the lowest-order term is proportional 
to n 7 /logn and dominates at time scales controlled by 
n <C h ~ exp(l/(l — 7)). Hence, for 7 ^ 1 and n 3> fi 
simple power laws are recovered while for 7=1 this log- 
arithmic term survives in the infinite time limit. The 
inset of Fig. 2 exemplifies this behavior explaining devi- 
ations between CTRW theory and the simulation results 
in the main figure. In other words, around z = 2 both 
normal and anomalous diffusion are suppressed due to 
logarithmic corrections in time leading to a vanishing K. 
We strongly suspect that such a behavior of the GDC 
is typical for dynamical phase transitions in anomalous 
dynamics altogether [6]. 

The stochastic CTRW theory presented above de- 
scribes only the coarse parameter dependence of the 
GDC. This motivates us to develop an approach that 
incorporates dynamical correlations. For this purpose 
we generalize the Taylor-Green-Kubo formula (TGK) 
[18] for maps to anomalous diffusion. We start by ex- 
pressing Eq. (2) via sums over the integer velocities 
v m = [x m +i\ — [x m ] [14]. However, in contrast to normal 
diffusion the anomalous dynamics generated by Eq. (1) 
is not stationary for z > 2 [6, 19]. This is reminiscent 
in the non-existence of an invariant probability density 
for this map, or mathematically, in the existence of infi- 
nite invariant measures [20]. Consequently, the resulting 
expression cannot be simplified using time-translational 
invariance, and we get 



K 



lim — 

n— too 77,7 



<H v l> +2< HY1 vkVk + i > 



k=0 



k=0 1=1 



(8) 

Focusing on the first term only, simulations confirm that 
it is proportional to n 1 . For z < 2 the dynamics is er- 
godic, hence this term is equal to n < Vq >, and we 
recover the random walk formula for normal diffusion 
[10, 17]. For z > 2 only generalized crgodic theorems 
may hold [20], and whether the CTRW result Eq. (7) 
can be derived from the first term is a non-trivial ques- 
tion. Considering this term as an approximation of K for 
all z, the numerical results are depicted in Figs. 1 (b) and 
2. The comparison of this data with the simulation val- 
ues based on Eq. (2) indeed shows that this term already 
provides a first step beyond the modified CTRW model: 
In the a dependence it reproduces the major irregulari- 
ties of K and even appears to follow the discontinuities 
discussed above. However, we cannot say yet whether it 
yields exact values for K at integer heights. 




FIG. 3: Comparison of the PDF obtained from simulations 
of the map Eq. (1) (oscillatory structure) with the analytical 
solution Eq. (10) of the fractional diffusion equation Eq. (9) 
(continuous line in the middle) for 2 — 3 and a = 8. The PDF 
was computed from 10 7 particles after n = 10 3 iterations. 
For the GDC in Eq. (10) the simulation result was used. The 
crosses (x) represent the numerical results coarse grained over 
unit intervals. The upper and the lower curves correspond to 
fits with a stretched exponential and a Gaussian distribution, 
respectively. The inset depicts the PDF for the map on the 
unit interval with periodic boundaries. 



Eq. (8) thus provides a suitable starting point for a sys- 
tematic understanding of the fractal GDC beyond CTRW 
theory. We remark that already the first term can be re- 
lated to de Rham-type fractal functions, which explains 
why it features irregularities [21]. Since the series ex- 
pansion in Eq. (8) is exact, working out further terms 
one must recover more and more structure in the GDC 
[10, 11, 14]. Note that the correlation function in the sec- 
ond term depends on two times, which enables to express 
K in terms of aged [8] de Rham-type fractal functions 
summing up to the exact value. 

We now turn to the PDFs generated by the map Eq. 
(1). The great success of CTRW theory derives from the 
fact that it correctly predicts both the power 7 and the 
form of the coarse grained PDF P(x, t) of displacements 
for a large class of models [5]. Correspondingly, the dif- 
fusion process generated by Eq. (1) is not described by 
an ordinary diffusion equation but by a fractional gener- 
alization of it. Starting from the CTRW model for the 
map Eq. (1) discussed above, one can derive the time- 
fractional diffusion equation 



d^P{x,t) _ 

dry = 



* / dP d 2 P 







or 



dx 2 



(9) 



describing the long-time limit of the PDF of the coarse 
grained dynamics with initial condition P(x,0) = 6(x) 
and D = KT(1 + 7 )/2, < 7 < 1. The fractional 
derivative is understood in the Caputo sense [22]. Time- 
fractional equations of such a form have already been 



extensively studied by mathematicians 



We remark 
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that the two other fractional diffusion equations pro- 
posed in Rcfs. [2, 8], which are based on a Riemann- 
Liouvillc fractional derivative, are equivalent to Eq. (9) 
under rather weak assumptions [24]. The solution of Eq. 
(9) expressed in terms of an M-function of Wright type 
[22, 23] reads 

p{x ' t} =^m; M H) (10) 

giving exactly the same asymptotics that was obtained 
in Ref. [5] for small and large values of £ = \x\j \f~Dt 1 1 2 . 
By using a series representation of M it can be demon- 
strated [24] that this form is equivalent to those expressed 
via H-functions [2] or one-sided extremal Levy stable dis- 
tributions [8], Fig. 3 demonstrates an excellent agree- 
ment between the analytical solution Eq. (10) and the 
PDF obtained from simulations for the map Eq. (1) if 
the PDF is coarse grained over unit intervals. However, 
it also shows that the coarse graining eliminates a peri- 
odic fine structure that is not captured by Eq. (10) [10]. 
This fine structure derives from the 'microscopic' PDF of 
an elementary cell (with periodic boundaries) as repre- 
sented in the inset of Fig. 3. The singularities are due to 
the marginal fixed points of the map, where particles are 
trapped for long times. Remarkably, that way the micro- 
scopic origin of the intermittent dynamics is reflected in 
the shape of the PDF on the whole real line: From Fig. 
3 it is seen that the oscillations in the PDF are bounded 
by two functions, the upper curve being of a stretched 
exponential type while the lower is Gaussian. These two 
envelopes correspond to the laminar and chaotic parts of 
the motion, respectively [25]. 

In conclusion, we have shown that the anomalous dy- 
namics generated by a paradigmatic subdiffusive one- 
dimensional map exhibits a fractal GDC under varia- 
tion of control parameters. The coarse dependence of 
this GDC and a non-trivial phase transition from nor- 
mal to anomalous diffusion are captured by a modified 
CTRW theory. Near the phase transition point the GDC 
is strongly suppressed by logarithmic corrections in time. 
A more detailed understanding of the GDC is provided 
by an anomalous TGK formula suggesting intimate rela- 
tions to aging, fractal functions and ergodic theory. The 
coarse-grained PDF of this anomalous dynamics is in ex- 
cellent agreement with the solution of a suitable frac- 
tional diffusion equation while on fine scales it reflects the 
microscopic details of the intermittent dynamics. Here 
we have only treated a subdiffusive map, however, we ex- 
pect these findings to be typical for spatially extended, 
low-dimensional, anomalous deterministic dynamical sys- 
tems altogether. Further studies will focus on a more 
detailed analysis of the anomalous TGK-formula, on a 
spectral analysis of this map and on treating superdiffu- 
sive dynamics along the same lines. 
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